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Abstract 

In this paper we analyze the transient behavior of the workload process in a Levy input queue. We 
are interested in the value of the workload process at a random epoch; this epoch is distributed as 
the sum of independent exponential random variables. We consider both cases of spectrally one¬ 
sided Levy input processes, for which we succeed in deriving explicit results. As an application we 
approximate the mean and the Laplace transform of the workload process after a deterministic time. 
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1 Introduction 

This paper studies the transient workload in a queue fed by a Levy input process X = {Xt}t>o', here 
the workload process, in the sequel denoted by {Qt}t>o, is defined as the reflection of X at zero. 
This workload process can be constructed from the input process X as the (unique) solution of the 
so-called Skorokhod problem, ll6l [T3lll4ll . It turns out that the process Q follows from X through 

Qt= Xt+ max{Qo, A}, 


where 


Ct := sup —Xs = — inf Xg- 

0<s<t 0<s<t 

The process {Ct\t>o is often referred to as local time (at zero) or regulator process Q. 
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As mentioned above, we are interested in the transient behavior of the workload process. In queue¬ 
ing theory transient analysis is a classical topic that is treated in various standard textbooks; see e.g. 
ll2ll5l lT2l . Typically, transient analysis is important in situations where the time horizon considered 
is relatively short, so that it carmot be ensured that the system is 'close to stationarity'. In addition, 
transient results are useful in cases fhaf fhe nef-inpuf process changes over time; if for insfance facil- 
ifafes fhe analysis of sysfems wifh fime-varying demand as well as fhe assessmenf of fhe impacf of 
specific workload confrol mechanisms. In general, fransienf analysis allows us fo assess fhe impacf 
of fhe initial sfafe Qq. 

The main contribution of fhis paper is fhe generalization of fhe existing resulfs on fhe fransient be¬ 
havior of fhe workload process {Qt}t>o- We consider n exponenfially disfribufed random variables 
Ti,..., T„ wifh paramefers qi,... ,qn and we analyze fhe joinf behavior of fhe vector 

(Qti, Qti+Tsi ■ ■ ■ j Qti + ...+T„) , 


with a specific focus on It is noted that this also directly yields Qt when T follows a 

Coxian distribution; see Section for some additional background on this claim. This observation is 
particularly useful owing to the fact that any distribution on the positive half line can be approxi¬ 
mated arbitrarily closely by a sequence of Coxian distributions (in the sense of convergence; see e.g. 
II 2 I Section III.4]). For the case of a spectrally positive input process, the results are given in terms of 
the Laplace-Stieltjes transform (LST), whereas we find an expression for the associated density for the 
spectrally negative case. 

Apart from the general results obtained, a second contribution lies in the reasoning behind our proofs. 
More specifically, our proofs reveal that the above formulas obey an elegant and simple tree structure. 
The transient workload behavior consists of 2” terms that can be recursively evaluated. We prove our 
results by induction; given that we know the expression for the quantity under consideration at n — 1 
exponential epochs, we derive the expression at n exponential epochs. In this induction step, from 
n — 1 to n that is, it can be seen how each term produces two offsprings, thus giving insight into 
the underlying structure. The idea behind the proofs yields a mechanism to address questions related 
to transient analysis at random epochs, which may help in obtaining a deeper understanding of the 
behavior of the underlying continuous-time queueing system. 

Transient analysis of queueing systems started with the analysis of the waiting times in the M/M/1 
queue ma. In lISlITSl the authors analyze the Laplace-Stieltjes transform of the waiting time process in 
the M/G/1 queue. The argument used there is also applied in ||6l, so as to derive Theorem |2.1 [ below 
for the case of a compound Poisson input process. The transient analysis of Levy driven queues is 
of a much more recent date; see e.g. ||2l|6l|8l for results on the workload process in a Levy-driven 
queue at an exponential epoch (which are briefly summarized in Section]^. As a direct application 
the authors of ID study clearing models, where special attention is paid to clearings at exponential 
epochs (relying on results on the workload at an exponential epoch in an M/G/1 setting). 
Concerning the structure of the paper, in Section 2 we present our notation, as well as the prelimi¬ 
naries that are needed in order to prove our results. In Section 3 we present the main results of the 
paper, which are Thms. 3.1 and 3.2 We support the final results with intuitive arguments based on 
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a tree structure; the proofs can also be interpreted along those lines. In Section 4 we present results 
obtained in numerical experiments. Section 5 contains the proofs of Thms. [3T| andIn all secfions, 
fhe specfrally posifive and specfrally negative cases are freafed separafely. Finally, Secfion|^confains 
conclusions and a brief discussion. 


2 Model, Notation, and Preliminaries 

In fhis section we presenf fhe workload af an exponential epoch for queues wifh specfrally posifive 
(Secfion |2.1| and specfrally negative (Secfion |2.2[ | Levy inpuf processes. These resulfs are heavily re¬ 
lied upon throughouf fhe paper, and in addition serve as a benchmark. In passing, we also introduce 
our notation. 


2.1 Spectrally positive Levy processes 

As mentioned, the building block of fhis paper is a Levy process X = {Xt}t>o- In case X is a spec¬ 
frally posifive process, henceforfh denofed by X e fhe Laplace exponent (j){a) := logEe““^^ is 
well defined for all a > 0. By applying Holder's inequalify we gef fhaf (/){■) is convex on [0, -|-oo) wifh 
slope ^'(0) = — EXi af fhe origin. In general, the inverse function '!/;(•) is not well defined and we 
work wifh fhe righf inverse 

ipiq) '■= sup{a > 0 : ^(o) = q}. 

For fhe case fhe driff of our driving process X is negative we observe fhaf = — EXi > 0 and 
fhus (/)(•) is increasing on [0, -l-oo). In fhis case fhe inverse funcfion is well defined. 

Our inferesf is in fhe transienf behavior of fhe workload process {Qt}t>o- We consider an exponenfially 
disfribufed random variable T wifh paramefer q (sampled independenfly from fhe Levy inpuf pro¬ 
cess) and focus on fhe fransform E^^ e““ where a > 0 and x denofes fhe initial workload. In fhis 
case fhe fransform E^; e““ is explicifly known, and is given in fhe following fheorem l^IHl ITSl . 

Theorem 2.1. Let X £ .5^+ and let T be exponentially distributed with parameter q, independently of X. For 
a > 0, X > 0, 

pOO 

E^ e-“ E^ 6"“ dt = 

Jo 

Using Laplace inversion fechniques HI, informafion abouf fhe process can fhen be inferred from fhe 
LST as if uniquely defermines fhe disfribufion of Qt, for each t and any initial workload x. 

2.2 Spectrally negative Levy processes 

For a spectrally negative Levy process X, henceforth denoted by X £ we define fhe cumulant 
4)(/3) := log E . This funcfion is well-defined and frnife for all /3>0, exacfly because fhere are no 
posifive jumps. We observe fhat $(•) has slope $'(0) = E Xi at the origin, thus 4’(/3) in general is no 
bijection on [0, -l-c»). We define fhe righf inverse fhrough 

'!'((;) := sup{/3 > 0 : 4>(/3) = q}. 


q - f{<y) 


f{q) 


=-■0(9)“ 
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When working with spectrally negative Levy processes, the so-called q-scale functions, ) and 

play a crucial role, particularly when studying the fluctuation properties of fhe reflecfed pro¬ 
cess mini. Forg > 0, lef lh(«)(a;), for a; > 0, be a sfricfly increasing and confinuous funcfion whose 
Laplace fransform satisfies 

/ = — -, /3>vl/(g); (2.1) 

Jo ^(P) - Q 

we lef equal 0 for a: < 0. From 0 Th. 8.1.(i)] if follows fhaf such a funcfion exisfs. Having 

defined fhe funcfion (•), we define fhe funcfion (•) as 

Z^i'lix) :=l + q f lh(«)(j/)dy. (2.2) 

Jo 

We immediafely see fhe imporfance of fhe g-scale funcfion in fhe densify of fhe workload process af an 
exponential epoch, given in fhe following fheorem ||6l Secfion 4.2], which is originally due fo Pisforius 

nni. 


Theorem 2.2. Let X G and let T be exponentially distributed with parameter g, independently of X. Lor 
a > 0, a; > 0 and /3 > 0, 


JQt G dy) = (^e '^^'^'>y'L{q)Z^'^'>{x) - qW^‘^\x - y)"j dy, 


and 


e-“ ®^da; =l( 


-F 


4'(g)-/3 a 


P \"^{q) + a <l>(/3) - g 4'(g) +aa + f 


The resulf on fhe LST of Qt in fhe above fheorem follows for /3 > 'F (g) by a direcf compufafion from 
fhe densify of Qt; by a sfandard analytic continuation argumenf fhe resulting expression fhen holds 
for any b > 0. 


3 Main Results 


In fhis secfion we presenf our main resulfs, viz. Thms. 3.1 and 3.2 In bofh subsections we firsf derive 
fhe workload behavior af two exponential epochs as fhis clearly demonsfrafes how fhe various ferms 
appear. Then we elaborafe on fhe mechanism for obfaining fhe workload af n exponential epochs, 
yielding an infuifively appealing free sfrucfure. The proofs can be inferprefed along fhose lines, and 
are given in full defail in Secfion]^ 


3.1 Spectrally positive case 

Suppose we have a spectrally positive Levy process X. We want to describe the behavior of the 
workload process {Qt}t>o at consecutive exponential epochs. We do this by considering exponen¬ 
tially distributed random variables Ti,..., T„ with distinct parameters gi,..., g„ and calculate, for 
Oi > 0 and some initial workload cc > 0, the joint Laplace transform given by 

E g —ai Cti —“2 Qti+T 2 +--- + Olri Qti + ...+T,i 
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It is instructive to first illustrate how to derive an expression for fhe joint transform af two exponen¬ 
tial epochs, i.e for E^, 2 t’i+t’ 2 . From Thm. 2.1 we have an expression for fhe fransform 

Consider now fwo exponentially disfribufed random variables wifh paramefers qi,q 2 - 


Then, conditioning on Qt^ in combinafion wifh applying Thm. 2.1 twice, yields 


POO 

E^ e"“i Ey 6"“" Vx{Qt, e dy) 

Jo 


92 


92 


92 - <^(a2) 
92 


92 - 4>icK2) 

E. - 




)) 


Ex{Qti G dy) 


9i 


^^( 92 ) 

ip p-(ai-|-y(g2)) Qti 

V^(92) " 

,-( 0 : 1 + 02 ) 2 ; _ Oil -|- 02 -tlj(qi)x 


92 - </>(a2) V91 - Hai + 02) V V'( 9 i) 

— _ 0^1 + V'(92) ^-i/i(i/i)a:\ \ p 2) 

'0(92) 9i - ^(ai + ^'(92)) V 0(91) // 

We see that by conditioning on the value of the workload at the first exponential epoch we can derive 
the transform at two exponential epochs. The above reasoning rests on the property that the process 
{Qt}t>o is a Markov process. 

Some special attention is needed for the case oi = 0 and qi = q 2 , i.e., when T has an Erlang-2 


distribution. From the last term in 13.21 we see that an additional limiting argument is required. 
A straightforward application of T'HopitaT then yields the expression for E^, e““ as in ftij Section 
4.1]. 


The main idea for the case of n exponentially distributed random variables Ti is very similar: condi¬ 
tion on the workload at the first exponential epoch, thus obtaining 


]g ^^1 “^2 Qti+T 2 —■■■—an Qti + ...+T„ _ 


' Ey e 


-as Cts eTs+...+T„ e dy). 

(3.3) 

to determine the transform at n exponential epochs 
given the joint transform at n — 1 epochs. At this point, it is useful to understand how the coefficients 
in the exponential terms of the transform appear; this is illustrated in Fig. [^below. Specifically, due 


Eqn. <|3.3[l is used in combination with Thm. 


2.1 


to the integration in < |3.3| and Thm. 2.1 it follows that each term produces two new terms when an 
exponential epoch is added (i.e., when moving from n — 1 to n exponential epochs), such that the 
transform at n exponential epochs consists of 2” exponential terms. We observe that in the expres¬ 
sion for n random variables the first term is, for every n, exp[—(oi a„)a:] (multiplied by some 

coefficient). The exponents exp[—(oi -I-... -I- q;;_i -I- tf}{qi))x], where I = 1,..., n, produce one expo¬ 
nential term of higher order exp[—(oi + ... + ai + tjj{qi+i))x] as well as one term corresponding to 
exp[—'0((3'i)x]; it is seen that the latter terms always appear at the 'even positions'. This mechanism 
is depicted in the tree diagram in Fig. where row n shows the 2” factors when we have n expo¬ 
nentially distributed random variables Ti,..., T„. For ease we only write the exponent at every node, 
hence the node oi -|-1/'(92) represents the term corresponding to exp[—(oi -I- '(p{q 2 ))x] (multiplied by 
some coefficient). In every row, the factors are counted from the left. 
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Qfl + 02 


^{qi 


ai + Ip{q2) 



ai+ip{q2) ipiqi) 


ai + a2 + 
^’( 93 ) 


^^(91 


i>{qi) ai+V'(92) V’(9i) 


Figure 1: The exponents in ED at every step 


We observe that the entire tree consists of subtrees starting from a node ipiqi) (aparf from fhe firsf 
element of every row). Suppose we have fhe elemenf exp[—(m + ... + ai-i + ip{qi))x] in the n-th 
row. This originates from a subfree generafed by an initial node tp{qi) thaf is / — 1 rows higher in fhe 
free. This follows from fhe facf fhat if we sfarf from fhe node ip{qi) we have fo move I — 1 times down 
and leff in order fo reach fhe node oi +.. . + a/_i+'!/’(9i)- So fhe node oi +.. . + ai-i+tj}{qi) in fhe n-fh 
row belongs fo a subfree spanned from fhe node il’iqi) iri the {n — I + l)-th row. For the ordering of 
terms, we assume that this initial node is at position 2j for some j = 1,2,..., 2"“*; we recall here that 
the nodes '(/’(gi) are located at the even positions of each row. Since the node is at position 2j there 
are 2j — 1 nodes in front of it. At every step downwards in the tree, the number of terms doubles 


since every term will give two new terms after using Thm. 2.1 Since we go down I — 1 rows, those 
2j — 1 nodes will produce in total {2j — 1)2^“^ = 2*j — 2*“^ nodes. Hence, we see that the element 
exp[—(oi + ... + ai-i + 'ilj{qi))x] in the n-th row is at the position 2*j — 2*“^ -1- 1. The numbering of 
the coefficients is based on this ordering. 



Figure 2: The coefficients of the exponential terms 

In Fig.|^ we use the following notation: 

o denotes the coefficient of the term exp[—(ni -|-... -I- a„)x]; 
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o 1 ) denote the coefficients of exp[—'!/'(qi)a;] (where j = 1 , 2 ,..., 2 " ^); 

o h| 2 'j- 2 '-i+i i) denofe fhe coefficienfs of exp[—(oi + ... + ai-i + ip{qi))x] (where I = 2 ,3,..., n 
andi = 1 , 2 ,..., 2 "-'). 


We nofe here fhat fhe superscripf (n) in fhese factors corresponds fo fhe number of exponential ran¬ 
dom variables considered (or, equivalently, in which row of fhe free we are). We now proceed fo fhe 
main resulf for fhe case of a specfrally posifive inpuf process. 


Theorem 3.1. Suppose we have n independent exponentially distributed random variables Ti, ..., with 
distinct parameters qi ,..., Qn- Then, for > 0 and x>t),we have 

n 

]g “^2 Qti+T 2 —■■■ — Otn QTi+T 2 + ...+Tn — I I _ ^.. .-\-OLji)x 

fJl Qi - Hoii + .. . + 


+ E E ( 3 . 4 ) 

1=1 j=i 


where the coefficients are defined below in Definition 


3.1 


The vecfors q = {qi,..., qn) and a = (oi,..., a„) are here explicifly included, so as fo show fhe 
dependence of fhe coefficienfs on fhe q's and a's. Lafer on fhese vecfors are omiffed fo keep fhe 
nofafion concise. 


Definition 3.1. For I — and j = 1,..., 2" '■we have 

n 

-(n) /- -\ TJ 


id al - TT _ "‘T _ TT __ 


where the Wi,n) given below in Lemma 

i = 1,2,... ,n — 1, are given through 


^(i,2h-2‘-Wl) ^ 


3.1 


^(n ,2 j 2 - 1 - 1 ) — 0 and the ^ for 


odd, 


even. 


r a,+i+d(*+h 2 h- 2 '-Wi) 

for 

2 i^._ 2 l-l_|_i 

2 * 

[ l/’fe-te) 

for 

' 2 h'- 2 '“'--|-l' 

2 * 


Remark 1. The terms are given from a recursive formula. The facf fhaf fhis recursion is well 
defined, follows because fhe lasf ferm equals zero (i.e., = 0 for all j's). 


Lemma 3.1. Consider j = 1, 2,..., 2" and take the binary representation of j — 1, i.e., j — 1 = 6o2° -I- 6i2^ -I- 
... -I- 6„_i2”“^. Then, for (or, equivalently, the sign of the j-th element in the n-th row of the tree 
presented above) we have 

g(i.n) _ (_2)Par{bo,. ..f>,x-l} 

where Par{6o, ..., dn-i} is 0 if the number ofl's in the binary expansion of j — 1 is even and 1 if it is odd. 


Similar fo fhe Erlang-2 sifuafion (i.e., n = 2 and <71 = qf), fhe case in which some of fhe qfs are fhe 
same has fo be freafed separately For insfance, n successive applicafions of 'I'Hopifal' lead fo an 
expression for E^; e““ , when T has an Erlang-n disfribufion. 
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3.2 Spectrally negative case 

In this subsection we concentrate on the case of a spectrally negative input process X. The joint 
workload density has a structure that is very similar to that observed for fhe LST in fhe specfrally 
positive case. Due fo the strong Markov property the joint density can be decomposed into 

IPa:(QTi G dj/i; • • • ; Qti + ...+T„ G ^Vn) = x{Qti G d?/i) • • • Py„_i(QT„ G dyn). 


That is, the joint density is simply the product of densifies af single exponential epochs, as given in 
Thm. 2.2 Henceforth, we focus on the density of the workload process at consecutive exponential 
epochs, i.e., 

Px(Qri+...+T„ G dy), y > 0. (3.5) 

First we illustrate how to obtain an expression for fhe densify Pa:(QTi+T 2 G dy) for some initial work¬ 
load X > 0 and y > 0 . From Thm. 2.2 we have an expression for the density G dy). Consider 

now two exponentially distributed random variables Ti,T 2 with distinct parameters qi,q 2 - Condi¬ 
tioning on Qti and applying Thm. 2.2 twice yields 

/*00 

Vx{Qti+T2 G dy) = / Pz(Qt2 G dy)Px(QTi G dz) 

J z^O 

pOO 


(-qiW^’i^\x -z) + P,(Qt2 G dy)dz 

nOO 

QiQ2 / — y)dz 

Jo 

pOO 

— / W^^^\x — z)Z^^^\z)dz 

Jo 

poo 

-q2'i/{qi)Z^‘^^'>{x) / e-'=^(«i)^lC(®)(z-y)dz 

Jo 

poo 

+ 'i>{qi)df{q2)e-'^''‘^^'>«Z^‘>^\x) / (z)dz 

Jo 


dy. 


(3.6) 


After some standard calculus and using the definition of the g-scale functions we find the expression 


':{Qti+T 2 G dy) — 


9192 (x-y)- d/{q2)qie-'^^‘^^^y ( x ) 

+ dy. (3.7) 

9i ~ 92 9i ~ 92 


Again the case qi — q 2 has to treated separately, by using I'HopitaTs rule; the detailed computations 
corresponding to this case can be found in Section 4.2]. 

We see fhat by condifioning on fhe value of fhe workload at the first exponential epoch we can derive 
the transform af two exponenfial epochs. As a nexf sfep our aim is to find an expression for | |3.5| for an 
arbitrary n > 0 and for exponenfially disfribufed random variables wifh paramefer qi{i = 1 ,..., n). 
Condifioning on fhe workload af fhe first n — 1 exponential epochs yields 

pOO 

Px(Qti+...+t„ G dy) = / P:j:(Qrj+.,.+T„_i G dz)¥^{QT„ G dy). (3.8) 

Jz^O 
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For the case of a spectrally positive input process (which was the topic of the previous subsection) one 
should condition on the value at the first exponential epoch, which allows the use of the induction 
hypothesis, but one needs to adjust the indices appropriately as the first exponential random variable 
is actually T 2 . For the spectrally negative case, however, conditioning on the value of Ti + ... + T„_i 
(and not only on Ti) allows us to circumvent this technicality. 

Moreover, the transition from step n — 1 to n can again be represented by using an elegant tree struc¬ 
ture that is similar to the one developed for the spectrally positive case. The expression for the density 
at n — 1 exponential epochs has 2”“^ terms and each term produces two new terms when integrated 
with the density (Qt„ G dy) (with respect to z). We also notice that in the expression for n exponen¬ 
tially distributed random variables the first term is always of the form * ... * ( 2 ; _ y) 

while the other terms are of the form * th(9i-i) * ... * (x), for I = 1,2,... ,n, multiplied 

by some coefficients that in general are functions of y. The underlying mechanism is illustrated in 
Fig. § In this tree the node — y) denotes the term ★ ... * FF(9 i)) ( 2 ; _ y'j while the 

nodes for I = 1,..., n, denote the terms * ff(9!-i) * ... * ( 2 ;). We see that at 

every row, say row k for ease, a new subtree with root * ... * ix) is created. 

These terms do not change as we move downwards in the tree since they only depend on the ini¬ 
tial workload x and do not take part in the integrations, similar to those carried out in | |3.6[ . Their 
coefficients change though, by a mechanism that is identified in the proof of our result. 


W^^^\x-y) 


Z(9i)(a 


W^‘i^\x-y) 








FF(®)(x-j/) Z(®)(a;) Z(®)(2;) Z^i^\x) Z(9i)(x) Z(9i)(a;) Z^i^^x) 

Figure 3: The convolution terms at every step 
We now proceed with the main result for the spectrally negative case. 

Theorem 3.2. Suppose we have n independent exponentially distributed random variables Ti,... ,Tn with 
distinct parameters qi,... ,qn- The density of Qti+...+t„/ given that Qq = x, is given by 


]Px(Qti+...+t„ G dy) = 


n 

(- 1 )” n 9 * • * ■ • ■ * ^^"'0 

n 2^~’- 

E E wi.o ^y) * ■ • ■ * 


1=1 j=i 


dy, 


where the coefficients t| 2 E 2 '-i+i are given in Definition 


3.2 
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Definition 3.2. For I — 1,... ,n and j = 1,..., 2” ^,we have the following expression 

n—1 


(v) - r^' 




1 


T ■ j Qi . / • 7\ I 1 Qm(j,l) 

— 1 , i—l i—m{j,l) + l ’ 


where m{j, 1) = min{fc € N : |fc—= !}• The terms ^ ^ +i.") are given below in Lemma 


3.2 


Lemma 3.2. Consider j = 1,2,..., 2” and take the binary representation o/2” — 2” — j = /3o • 2° + ... + 

Pn-i ■ 2"“^. Then, for fT”.) equivalently, the sign of the j-th element in the n-th row of the tree presented 
above) we have the following formula 

(fTn) _ ^_2^^Par{/3o,/3i,...,/3„_i} 

where Par{/3o, ■ ■ ■, /3n-i} is 0 if the number ofl's in the binary expansion o/2" — j is even and 1 if it is odd. 


Using the result obtained in Thm.|3^we can find an expression for fhe fransform wifh respecf fo fhe 
initial workload as well; again analytic continuation is used fo obfain fhe resulf for any /? > 0 . 


Corollary 3.1. For a > 0, f3 > 0 and for n independent exponentially distributed random variables Ti,..., T„ 
with distinct parameters qi,... ,qn, we have 




n 


+ P ^(/3) - qi 


1=1 j=i 


(2‘j-2‘~'^ + l,n) 




n 


n—1 

n 


i—l, 


■ I qi qi+1 ■ , ■ , ,qm(j,i) 

*=i z=m(j,i) + l ' 


qi+1 


a + Ti{qm(j,i)) H $(/3) -q^ p 


where m{j, 1) and c 


,(2'j-2'-i + l.n) 


are given in Definition 3.2 and Lemma 


3.2 


The case in which some of fhe qfs are fhe same should be freafed separafely again. For insfance, 
fhe densify of VxiQr € dy) for T having an Erlang-n disfribufion follows affer n applications of 
I'HopifaTs rule. 


4 Numerical Calculations 

In fhis section we presenf numerical illusfrafions of fhe fransienf workload behavior. We consider 
examples corresponding fo fhe specfrally positive case (noting fhaf fhe specfrally negafive case can 
be dealf wifh similarly). The expression found in Thm. [3T| is, from an algorifhmic sfandpoinf, highly 
attractive; fhe only drawback is fhaf for every n we have fo compufe 2 " ferms, fhus increasing fhe 
compufafion fime significanfly af every sfep. In our illusfrafions, we consider fhe impacf of n, i.e., fhe 
number of exponenfial variables. We also commenf on ways fo defermine fhe workload disfribufion 
af a fixed (deferminisfic, fhaf is) fime; fhe mean idea fhere, as we poinf ouf in more defail below, is fo 
approximafe a deferminisfic epoch t by fhe sum of exponentially disfribufed random variables wifh 
appropriafely chosen paramefers. 
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We focus on two specific Levy processes: Brownian motion and the Gamma process. For the case 
of Brownian motion, the input process X is a Brownian motion with a drift, henceforth denoted by 
X G a^). Then, (j){a) := logE = —ad + q:^(t^/ 2, and the right inverse function is 


-iPiq) = 


d + + 2cr^q 


For reflected Brownian motion (i.e., the workload of a queue with Brownian motion as input) there 
is an explicit expression for the conditional distribution P(Qt < y\Qo = x), for y > 0, see e.g. El 
Section 1.6]. It is a matter of straightforward calculus to use this formula to find an expression for the 
transform 

E, e-“ 2* = / e-“^ P(Qt e dy| Qo = ^). 

do 

This result is used to evaluate the performance of our procedure in case a fixed time t is approximated 
by the sum of exponentials. 

The Gamma process is characterized by the Levy-Khintchine triplet (d, cr^, 11), where cr^ = 0; the Levy 
measure is given, for some /3 ,7 > 0 , by n(da;) = {fl/x) , for a; > 0 , and the drift is d = xn(da:). 
From the definition of the Levy measure we see that the Gamma process is a spectrally positive 
process with a.s. non-decreasing sample paths. We also add a negative drift such that the Laplace 
transform is equal to 


(p(a) := logEe = /31og 


a 


pa. 


where p > /S/y in case of a negative drift d = {P/'y) — p. For the Gamma process with parameters 
7 , /3 > 0 and a drift (/l/y) — p we use the notation G( 7 , /3, p). If the input is a Gamma process there is 


no explicit expression for the transform E;j 


o-a Qt 


in contrast with the case of a Brownian input. 


Suppose now that we wish to characterize the distribution of Qt for a deterministic t. The idea is 
that we can approximate t by a sum of, say n, independent exponential random variables. An op¬ 
timal choice of the parameters qt then follows from solving the following constrained optimization 
problem: 

n ^ n 

min Var(Ti-I-...-I-T„) = min^^ “2 s.t. ET^ = t. 

i—l 2—1 

This constrained optimization problem has solution qi = ... = qn = n/t. A complicating factor is 


that in Thm. 3.1 the parameters qt should be chosen distinct. To remedy this, we propose to impose a 


small perturbation of the optimal g/s such that they are distinct: 

a,), (4.1) 

where the at are suitably chosen small numbers that sum up to 0. In the two tables that follow we 


- = -(l 

Qi n 


present the numerical results obtained from calculating the expression in Thm. 3.1 for the case X G 
Bm(—1,1) (Table0 and for the case X G G(l, 1,2) (Tablej^. Here, we consider the situation of a: = 0 
and t = 1. The parameters qt are chosen according to <|4.1| with, if n is even, the a/s given by 


Qj — 


0.01 G iff =!,...,§ 
— 0.01 ■ i if i = § -I- 1 ,..., n. 
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Table 1: Numerical approximations for X G Bm(—1,1), a: = 0 and t = 1 



n = 1 

n = 4 

n = 6 

n = 7 

n = 8 

exacf value 

relafive error 

a = 0.1 

0,9647 

0,96064 

0,96021 

0,96005 

0,96001 

0,95914 

-0,09 % 

a = 0.2 

0,9318 

0,92410 

0,92327 

0,92299 

0,92300 

0,92128 

-0,19 % 

a = 0.3 

0,9011 

0,89008 

0,88892 

0,88851 

0,88836 

0,88611 

-0,25 % 

a = 0.4 

0,8723 

0,85836 

0,85688 

0,85638 

0,85608 

0,85338 

-0,32 % 

a = 0.5 

0,8453 

0,82870 

0,82696 

0,82637 

0,82590 

0,82285 

-0,37 % 

a = 0.6 

0,8199 

0,80094 

0,79896 

0,79828 

0,79786 

0,79432 

-0,44 % 

a = 0.7 

0,7960 

0,77488 

0,77270 

0,77196 

0,77237 

0,76760 

-0,62 % 

a = 0.8 

0,7735 

0,75040 

0,74803 

0,74723 

0,74625 

0,74254 

-0,50 % 

a = 0.9 

0,7522 

0,72735 

0,72482 

0,72397 

0,72415 

0,71900 

-0,71 % 

a = 1 

0,7321 

0,70562 

0,70295 

0,70205 

0,70205 

0,69684 

-0,74 % 


(If n is odd we choose a„+i /2 = 0 and fhe resf as indicafed above.) In fhe firsf fable, X G ]Bm(—1,1), 
we fake n = 1,4, 6 , 7 ,8 and compare our approximations with the exact values obtained from E^; e““ 
for differenf values of a. In fhe lasf column we presenf fhe relafive errors befween fhe exacf value and 
fhe approximafion value for n = 8 . 

If should be realized fhaf fhe numerical procedure has ifs limifafions. Firsf, from fhe expression in 
Thm. 3.1 we see fhaf af every sfep we have fo compufe 2" ferms, which complicafes fhe compufafion 
for n large. We also see fhaf when fhe paramefers are 'almosf equal' i.e., in < |4.1| is small) we 
add and subfracf ferms fhaf are large in absolufe value (as fhe denominafors feafuring in fhe resulf 
of Thm. 3.1 are close fo zero), which pofenfially causes insfabilify. Our numerical fesfs show fhaf fhe 
choice of fhe paramefers qi influence fhe numerical sfabilify; for the parameters indicated in | |4.1| the 
results begin to deviate for n > 9 due fo numerical issues. From Table corresponding fo fhe case 
of a Brownian inpuf process (for which we can compare wifh exacf resulfs), we see fhaf for n = 8 
our relafive error is below 1%. For fhe case of a Gamma inpuf process, we verified fhaf fhe fransform 
converges fo fhe sfeady-sfafe workload as given by fhe generalized Pollaczek-Khintchine formula fH 
Thm. 3.2]. 

As a second application of Thm. 3.1 we use fhe resulfs obfained in Tables and fo approximafe 
fhe value of E^^ Qt for fhe two cases X G Bm(—1,1) and X G <G(1,1,2), essentially relying on nu¬ 


merical differentiation. By considering an a sufficiently small and an n sufficiently large we use the 
approximation 

Ea: Qt 


1 - Er - 


a 

We present our findings for fhe cases X G Bm(—1,1) and X G G(l, 1,2), respectively, displaying fhe 
qualifafive behavior of E^; Qt as a funcfion of fime for various values of x. For fhe mean value of 
fhe sfafionary workload we know fhaf E Q = as follows direcfly from fhe generalized 

Pollaczek-Khinfchine formula. 
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Mean value 


Table 2: Numerical approximations for X G G(l, 1,2), a; = 0 and t = 1 



n = 1 

n = 4 

n = 5 

n = 6 

n = 7 

n = 8 

a = 0.1 

0,97582 

0,99046 

0,99037 

0,99032 

0,99028 

0,99026 

a = 0.2 

0,95527 

0,98148 

0,98130 

0,98121 

0,98112 

0,98108 

a = 0.3 

0,93754 

0,97300 

0,97275 

0,97261 

0,97249 

0,97243 

a = 0.4 

0,92205 

0,96499 

0,96465 

0,96448 

0,96432 

0,96425 

a = 0.5 

0,90838 

0,95739 

0,95699 

0,95678 

0,95659 

0,95662 

a = 0.6 

0,89621 

0,95018 

0,94972 

0,94948 

0,94925 

0,94936 

a = 0.7 

0,88530 

0,94333 

0,94281 

0,94254 

0,94228 

0,94201 

a = 0.8 

0,87543 

0,93681 

0,93623 

0,93593 

0,93565 

0,93565 

a = 0.9 

0,86647 

0,93060 

0,92996 

0,92964 

0,92933 

0,92885 

a = 1 

0,85828 

0,92467 

0,92398 

0,92363 

0,92330 

0,92273 


Mean value approximation for Brownian motion 



Figure 4: Mean value approximation with n = 7 
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Figure 5: Mean value for Gamma process 

In Figs. 1^ and we observe three different scenarios corresponding to different values of the ini¬ 
tial workload. When the initial workload is 0, the mean workload increases and converges to the 
mean value of the steady-state workload. This follows directly, as, for any Levy input process, 
Qt = supo<s<j Xs when the initial workload is 0, implying that Eq Qt is increasing over time. When 
the initial workload is slightly above the steady-state workload, it is interesting to notice that E Qt 
first decreases below the steady-state version, and then converges from below. For higher initial 
workloads, E Qt is always decreasing and converges to the steady-state value from above. 


5 Proofs 

In this section we prove the main results. For both cases of spectrally one-sided processes, we first 
show an auxiliary lemma relating to the signs of each term. The main results are then proved using 
induction. 


5.1 Proof of Theorem 13.11 


Before deriving the main result, we first prove Lemma 3.1 which gives the sequence of the 2" signs 
that appear in the expression of the transform at a time epoch corresponding to the sum of n exponen¬ 
tially distributed random variables. From Thm. 2.1 we see that for n — 1 the signs of the coefficients 
are -I-, — . For n = 2 and from Eqn. 13.21 we see that the signs are -I-, —, —, -I- (where it is noted that we 
use the ordering of the terms presented in Fig.|^. Since we know how the terms are produced when 
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we go from the step with n exponential times to the step with n + 1 exponential random variables 
(see Section [3d) we see that the signs at every step can be represented again by a tree graph. In this 
tree, row n again consists of 2" nodes and, starting from the left, the nodes represent the sign of every 
factor when the expression is written as in Eqn. (|3.4). 



+ 


/\ 


+ - 





/\ 


- + 


Figure 6: The sequence of the signs at every step 


We see that row n + 1 can be derived from row n when substituting every + in row n by the pair 
+, —, and every — by the pair —, We can understand why this holds by looking at the expression 

Denote by the sign of the j-th element 


2.1 


and the mechanism analyzed in Section 


3.1 


in Thm. 

in the n-th row in the above tree. Then for j = 1,2,..., 2", corresponds to the sign of the j-th 
coefficient when considering n exponentially distributed in Eqn. ||3.4). 


Remark 2. We observe that, because of symmetry, for the signs of the fc-th row it holds that, for 

^U,k) ^ 


Hence the signs j and j' + 2^ ^ in the fc-th row will always be opposite. 

Proof of Lemma |3d| We prove the lemma by induction on the number of exponentially distributed 
random variables. 


(i) For n = 1 we have two nodes and this case corresponds to the signs of the expression derived 
for one exponentially distributed random variable T. We have that = +1 


in Thm. 


2.1 


and = —1. Then we need the binary expansions of 0 and 1 which have no I's and one 1, 
respectively. We see that = (—1)° = 1 and = (—1)^ = —1. 

(ii) We assume that the lemma holds for n = k. Hence, for j = 1,..., 2^, we have 


Here we make the following observation. In the tree presented above, consider an arbitrary row 
n. The 2” signs of that row and the first 2" signs of the (n -I- l)-th row are the same. 

Now consider the (fc-l- l)-th row. Using the observation above and the induction h 5 rpothesis, the 
lemma holds for the first 2^ signs of this (fc + l)-th row. Hence, we need to prove this statement 
only for j = 2^ -I- 1,..., 2^+^. 
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For j = 1,2,..., 2^ we have 

j,(i.fc+l) _ ^_ 2 ^Par{ 6 o,...,bfc_i} _ ^^Par{f)o , 0 } 

where j — l = bQ + bi-2 + ... + b^-i ■ 2^“^ + 0 • 2^. Consider now the element j' = j + 2^. 
From Remark|^we know that c*--' ’^+i) = —ch.^^+i), vVe also know that the binary expansion of 
f has one more 1 than the binary expansion of j since we add 2 ^, i.e., / — 1 = j — 1 + 2 ^ = 
bo + .. . + bk-i - 2 ^“^+ 1 - 2 ^, which shows fhat (_i)Par{bo,--,b)=-i.i} = _^_]^^Par{bo,...,bfc_i} leading 
fo 

_ ^_ 2 ^^Par{bo,...,bfc} 

for all j = 1 , 2 ,..., 2 ^+^. 


□ 


Before proceeding wifh fhe proof of Theorem 3.1 we presenf some general remarks which are used in 


fhe proofs of Thms. 3.1 and 3.2 


Remark 3. For I = 2,... ,n and j = 1,, 2"“* we observe fhe following 

(a) 2^j — 2*“^ + 1 is an odd number. 

(b) For alH = 1, 2,..., 1 — 2, 


2'j - 2 ‘-^ + 1 
2i 


2'-i-(2j_l) + l 


= 2 ‘-*-'( 2 j-l) + l, 


which is always an odd number. In addifion, 

r 2 *j _ + 1 


= 2j 


is an even number. 


(c) 


(d) For i = 0,1,... 


2^j - 2} 


2 '-i 11 


r 1 11 

' 2'' 

— 

^ ~ 2^ ¥ 


= J 


■2'j - 2'-i -6 1' 


'3 

1 1 


'3 ' 

2l+i 


2* 

2^+1 ' 2 *+® 


2* 


Proof of Theorem [3T| We use induction on fhe number of exponential random variables Ti,..., r„. For 


the proof if is sufficienf fo sfarf wifh n = 1 (where if can be readily checked fhaf Thm. 3.1 holds for 


n = 1), buf fhe case n = 2 is more insfrucfive. The joinf fransform for n = 2 can be found in Eqn. 13.21. 
Firsf of all, when n = 2 we have in tofal 2^ = 4 terms. We see that the even terms correspond to 
exp[—'ijj{qi)x], and the third term corresponds to exp[—(oi + 'ip{q 2 ))x]. According to ( |3.4| the coeffi¬ 
cient of exp[—(oi -I- 02 ) 2 ^] musf be equal fo 


12 


qi 


92 - 0(Q!2) 9i - fiai + 02 ) ’ 
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following directly from < |3.2| . We have two coefficients corresponding to exp[—which according 
to (3.4 1 should be equal to ^[ 2 ^) Using Definition 3.1 we find the following expressions 


"S.-n— 


n 




92 


9i 


Ol + 02 


g,2 - ^f)(a2) 9i - + 02 ) ^{qi) 

as dU’ 2 ) = ip{qi), = 02 , and = 0. Moreover, 


U’.,=n- 


9i 


Qi — (/)(ai + 11 dii-iA) 


n 


- d(*d) 


where we see from the table for the factors (500 Definition 3.1 1 that rfUd) _ = ip{q 2 ), 

and dUd) = Q. This leads to the following result 




92 


02 


9i 


oi + ' 0 ( 92 ) 


dd) 52 _ 0(q!2) ■ 0 ( 92 ) 9i - 0(ai + 0 ( 92 )) 0 ( 91 ) 

For the last term, the coefficient of e“di+V’(92))a:^^ y^Q g 0 |- 


= - 


9i 


92 


Ud) - (j)(^ai + 0 ( 92 )) 92 - 0 ( 02 ) dh h3) 


n 




dh.3) 


Since d^d) = 0 ( 52 ) and d^d) = 0 , this agrees with Eqn. (3.41, and thus the results holds for n = 2. 
We now assume that our formula holds for n = k — 1. Hence we have that 

k-l 




q^ 


^-(ai+...+ak-i)x 


qi - 4>{ai + ... + Ofc-i) 


k-l 2 '=-'-! 


r(fc-l) 


^-(ai + ...+ai_i+y(qi))x 


(5.1) 


1=1 i=i 


where the coefficients U( 2 q- 2 '-i+i i) given by Definition 3.1 for n = fc — 1 and the signs of all the 
factors are given by Lemma 3.1 In the induction step we prove this theorem for n = k given that it 
holds for n = k — 1. The expression for n = k is derived from calculating the integral 

poo 

C := e"“i -■■■-“'= 2 ti+...+t;, = / g-a 2 Qt^--- o^k Qt^+.-.+t^ VxiQxi G dy), 

Jo 

where the expectation in the integral is known by the induction hypothesis. Here we see that we 
must raise all indices in •ED by one when we do the calculations because we start from time T 2 with 
parameter q 2 instead of from Ti. Combining the above with (|5.1[l, we obtain 


fc-i 


^ = n 


9i+l 


9i+l ~ 0(Q^i+l + + Oik) 


^-{ai + ...+ak)y 1 


c(Qti G dy) 


k—i 2 '=-!-* 

r(fc-l) 

Z^ Z^ -^(2h-2'-i + l,/) 

1=1 j=l 


^-{ai+...+ai+tp(qi + i))y i 


Gd 2 /)=:Z + IT (5.2) 
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The two integrals in < |5.2| can be computed using Thm. 2.1 Each integral gives two new terms, cor¬ 
responding to a move down and left for fhe firsf ferm, and down and right for the second term in 
the trees presented in Figs. and The exponents are easily observed after an application of Thm. 


2.1 Therefore, below we primarily focus on fhe coefficienfs. When considering such infegrals fhe 


two ferms obfained are referred fo as fhe firsf and second ferm and are denofed by adding a 1 or 2 as 
indices fo I and XX. We now successively consider (fhe coefficienfs of) X^, XX^, X^, and XX^. 


o Coefficient ofXi. The coefficienf of exp[—(oi -I-... -I- ak)x] is found, using Thm. 2.1 from fhe firsf 
ferm of fhe infegral 

/•OO k — 1 


r,-otiy 




/O 


n - (j){a^+i + ... -I- Ofc)' 


^-ia2 + ...+ak)v 1 


c(Qti e dy), 


which is 


% 


<7i 




n + ... -t Ofc) qi - (j){ai + ... -|- Ofc) Qi - ficti + ... -t Ofe) ’ 


n 


This corresponds fo fhe coefficienf of fhe firsf ferm in Thm. 3.1 


(5.3) 


can 


o Coefficient ofXXi. For Z = 2,3,..., fc if is seen fhaf fhe ferms for j = 1,2,..., 

be derived from fhe ferms i_i) ^>7 taking the first term of fhe infegrals (fhis corresponds 

fo a move down and leff when we look af fhe free in Fig.|^: 


L 


(fe-i) 




(5.4) 


From Thm. I2.1l we obfain 


-(fc-i) 


qi 


T i^) _ r f _ 


= ^( 2 ‘-ij- 2 ‘-'’ + l.fe-l) 


k-1 


Qi+1 


k-l 

n 


Oii+l 


ryl-l - rfl-2 




qi 


i—l — 1 

r)Z-l • r,l-2 


d(i,2'-ii-2'-2 + l) 
k 


,.( 2 ‘-ij- 2 ‘-'' + l.fc-l) , 


n 


qi - 4>{ai -I- ... -t ai-i -t f^qi)) 
qt 


1^2 
r\l — 1' r\l — 2| 




a, 




qi 


n dh-1.2' ij-2' 2_i_i) _i_ ’ 

where j = 1,2 ,..., 2'^“^; here Jd. 2 ' ^j- 2 ' ^+i) jg given by 

is odd, 

. ol-l • r,l-2 


J(i, 2 '-ij- 2‘-2 + l) ^ 


r a,+i + 

if 

2 i-lj_2 i-2_|_i 

2i-l 

[ V'fe-hi) 

if 

'2'“ij-2'“^ + l' 

2i-l 


IS even. 
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This table follows from Definif ion |3 .1 [ and fhe observation fhaf fhe facfor c?*^*’^* initially was 

fhe facfor added fo fhe ferm (fhis is why we use fhe nofafion d for fhese ferms); fhis is due fo fhe 


facf fhaf in 15.41 all indices are raised by one. In order fo bring fhis info fhe form of Definifion 3.1 
observe fhe following: 


we 


(a) Concerning fhe signs we have fhe relation +i.fc-i) _ ^(2 j -2 +i,fe) £qj. i _ 2 ^ _ ^ ^ 

and j = 1,..., We see fhis as follows. From Lemma 3.1 we see if is sufficienf fo show fhaf 

fhe numbers — 2*“^ and 2*j — 2^“^ have fhe same parify. Buf fhis holds as 2*j — 2*“^ = 

2{2^~^j — 2*“^). Infuifively we can see fhis from fhe free graph in Fig.|^ every fime we move 
down and leff fhe sign is always fhe same. 

(b) Concerning fhe labeling of fhe ferms, using fhe facf fhaf 


■2'-ij - 2'-2 + 1' 


'2^j - 2'-i + 1' 

to 

T 


2* 


(which we obfain from Remark]^ we obfain 

^p,2h-2'-i+l) ^ J(i,2'-L-2*-Wl)^ 


(5.5) 


(c) From fhe four properfies in Remark|^we see fhaf 




l,2‘j-2‘-^ + l) _ 


— 02 + . . . + 0;_1 + Tpidl) ^rid 


= 1 , 


^(i-l, 2 h- 2 '-i + l) 

for alH = 1, 2 ,..., Z — 1. 

The argumenfs in (a)-(c) show fhaf, for I = 2,3,... ,n, j = 1,2,..., 2^“*, 


T (^) 




n 


di 


n 




h2*i-2'-i + l.Z) - 11 _^^P.2h-2‘-i + l)) 11 rfti-1.2h-2'-i + l) 

i—1 ^ ^ '' i—i 

where fhe are given by fhe fable in Definifion; 


3.1 


o Coefficient 0 /X 2 . For fhe ferms j = 1,2,..., 2^~^ (i.e., fhe coefficienfs of e^p[—'ip{qi)x] for 

k exponentially disfribufed random variables) we observe fhaf fhese are given from all ferms in fhe 
previous sfep, one from each (fhis corresponds fo moving down and righf in fhe free graph in Fig. 
or Fig.j^. The firsf ferm, resulfs from fhe infegrafion 


di+l 


, k — 1 

^ Qi+l ~ f{o^i+l + . . . + Ofc) 


e-(“^+-+“'')^P.(QTi edy). 


which leads fo 


4:1)=-n 


qi 


ai 


O-k 


Qi - <p{ai + ... + ak) qi + f{ai + ... + Ofc) -ip(qi) 


Since I = 1 and j = 1, we have for i = 1,2,..., k 

2 

¥ 


= L 
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showing that dh, 2 ) _ as- Furthermore, we see that for alH = 2,3,..., fc 


ai + (fh’i) 
dh-1.2) 


and, hence, we get 

Oi + _ oi + _ oi + ... + Ofc 

IJ- dh-1.2) “ d(0,2) “ ■ 

By using these facts, it follows fhaf 

j-(k) __TT^__ TT Qj + 

~ + dh’2)) dh“i’2) 


corresponding fo Definifion 3.1 and Lemma 3.1 


(5.6) 


o Coefficient of IJ 2 - In general, fheferms L| 2 iVij_ 2<+2 i)'^ ^ ^ 

are derived from fhe infegrals 



T (^“ 1 ) 


g-«iy 


e 


(a2 + ...+a,+V-(9, + i))yp^(g^^ ^ 


(5.7) 


Consider fhe ferms h|2!+ij_2i+2 i) for 1 = 1,..., A: — 1 and j = 1,..., 2^ ^ ^ 
weobfain,for I = I,... ,k — 2 and j = 1, 2,..., 2^~‘^~^, fhaf 


From fhe infegral in (5.7 


T (^) 

(2'+ij-2'+2,l) 


k-1 

_ (2h-2'-^+l,fe-l) TT g»+l 

■ l{ %+l - '/'(Oz+l + d-h+h2h-2'- + l)) 

W Oj+i + _gi_ Oi + ... + ai + V>(g;+i) 

11 J(i.2'i-2'-i+i) _ ^(q,^ Q,^ _|_ 


(2'j-2'-i + l.fe-l) TT _ % _ 

iigz-0(a* + d'h.2h-2‘-+i)) 

A Qfj + d(^’ 2 h- 2 ' ^+ 1 ) _ ^ _ Oi + . . . + g; + fjqi + l) 

11 J(*-i, 2 'i- 2 '-i+i) + ... + ai +'f{qi+i)) f{qi) 


where fhe facfors d6’2 i 2 +1) ^j-g gjyen by 

J(j.2h-2'-i + l) ^ 


r g,+i+d'6+i’2h-2'-Ui) 

if 

'2h'-2'“'^ + l' 

2i-l 

[ V'(ft+i) 

if 

’2'j-2'“'^ + l 

2i-l 


is odd. 


IS even. 


Using fhe same observation as in (5.5 <, it is found, for j = 1,..., 2^ ‘ ^ and i = I + 1 ,..., k, fhaf 

j(i,2' + U-2'+2) ^ J(i,2h-2'-i + l)^ 


From Remark |^(a)-(c) we see fhaf 

^(1.2‘+U-2'+2) = + , . . + rf(0,2'+U-2U2) ^ ^(q^) 
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and, for * = 2, 3,..., Z, 


^(i-l,2* + ij-2'+2) 


= 1 . 


These observations allow us to write ^|2'+ii-2'+2 i) follows 


-(fc) 


^(2i+ij-2'+2.1) 


= _„(2h-2'“^ + l.fe-l) . 


n 




l^ig^_^(a^+(ih,2‘+ij-2'+2)) ^(i-l,2< + ij-2‘+2) ■ 


n 


a 


^+rfM' + ij-2'+2) 


Concerning the signs, we obtain the relation ^ ^ +2,fe) = _£;(2 j 2 +i,fc 1) gjj-((,g ^]-(g numbers 
2*+^ j — 2^ + 1 = 2(2* j — 2*“^) + 1 and 2*j — 2*“^ have opposite parities. This final expression agrees 


with those presented in Thm. 3.1 


Now, we combine the above results to complete the proof. Using the coefficients of Ii and I 2 , i.e., 
x3) and (|5.6b, we can rewrite ||5.2b to 


^=n 


Qi 


f=i ~ ■■■ 


^-{ai + ...+ak)x _ 




Ui + 

qi-4>{.ai + 

2=1 ^ ^ ' 2 = 1 


n 


n 


-ip(.qi)x 


n cyK — l 

^ _ pOO 

+ E E ^&-E2 '-Ni,;-i) / e dy). 

7 n ^ — 1 ^ 0 


Z=2 i=l 


ffe) 

Using the definition of T(2 p/ in conjunction with the coefficients IIi and II 2 and Definition 
above expression can be written as 


3.1 


the 


^=n 


__ (^a^ + ,,,+ak)x _ r{k) 

fji ■■■ 


k 2^- 


k 2^ 


I p-(«i + ...+«!-i+'0(9i))a; I r(fc) p-i>(qi)x 


/=2 i=l 


i=2 i=l 


It remains to write the last sum in the desired form. This double sum has in total 2*““*^ — 1 terms, 
and we observe that for I = 2,... ,k and j = 1 ,..., 2*^“*, 2^j — 2*“*^ + 2 defines a partition of the even 
numbers 4,6 ,..., 2* irito k — 1 classes each one containing 2*^“* numbers. Relabeling the terms with 
only one subscript, we can write this double sum as 


k 2“ 


E E '^(2h-2'-i+2,l)' 


-‘>p{qi)x 


f 

;=2 2=1 i=2 


= E 


(k) -i){qi)x 


1)' 


where i = 2 j — 2 ^ + 2 for I = 2,..., k and j = 1,..., 2*^ ■ From the above it follows that C can be 


written as the expression in Thm. 3.1 This completes the proof. 


□ 


5.2 Proof of Theorem 13.21 

For the spectrally negative case we use a tree graph to illustrate, similar to Fig.|^ how the coefficients 
of the convolution terms change from step n to n + 1. Based on the reasoning that led us to the 
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tree diagram in Fig. for the spectrally negative case, we adopt the numbering of the coefficients 
asinFig.g 


L 


( 2 ) 

1 






( 1 ) 






L 


( 2 ) 

4,1 



Figure 7: The coefficients of the convolution terms 


Comparing this tree graph with the one corresponding to the spectrally positive case (i.e., Fig.|^, we 
observe that now the numbering of fhe coefficients is 'mixed'. In the spectrally positive case every 
time we went right a new subtree with root V'('Zi) was generated, whereas here only the first time 
we turn right a subtree with root * FFti'-d * ... * (^x) is generated. So, in the spectrally 

negative case, convolution terms of the same order gather in one subtree. When moving downwards 
in the subtree the coefficients change according to a recursive pattern. This mixed enumeration of 
fhe coefficienfs is due fo fhese fwo characterisfics. In the spectrally positive case terms of all possible 
orders are generafed in all subfrees generafed by some ip{qi). From Fig.j^we see fhaf the coefficient 
■^( 2 'j- 2 '-i+i.i) generates the coefficients ^( 2 *j- 2 *-i+i+ 2 ",i)' observe 

that, for 1 = 1, ..., n, we have a fofal of ferms of type * ... * (x), which is 

why we choose to label the coefficients of these terms by the numbers — 2*“^ + 1. This trick leads to 
an expression which is relatively easy to work with in the proof, and af fhe same time has a structure 
similar to the one featuring in Thm. 3.1 for fhe spectrally positive case. 

Let us now first consider the sign of fhe j-fh term when we have n exponentially distributed random 
variables. From Thm. 2.2 we see that, for n = 1, the signs of fhe coefficients are —, +. For n = 2 and 
the expression in Eqn. < |3.7| it turns out that the signs are Following how the terms are 

produced when we go from fhe sfep wifh n exponential times to the step with n +1 exponential times 
(Section |3.2| and Fig.j^, the signs at every step can be represented by the tree graph in Fig.|^ In every 
row, starting from leff to right, the nodes represent the sign of every factor when our expression is 
wriffen in fhe form presented in Fig.|^ 
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/\ 


+ - 





Figure 8: The sequence of the signs at every step 


We see that row n + 1 can be obtained from row n if we substitute every + by the pair —, +, and every 
— by the pair +, —. This holds due to the mechanism analyzed in Section 3.2 i.e., the order in which 
the integrations are carried out. Denote by the sign of the j-th element in the n-th row in the 
tree. Then j = 1,2,..., 2" and corresponds to the sign of the j-th coefficient when we have n 


exponentially distributed random variables in the expression considered in Thm. 3.2 


Proof of Lemma [T2| We prove this lemma by induction. 

(i) For n = 1 we have to find the values of and For j = 1 we need the binary expansion 
of 2^ — 1 = 1, which has one 1, while for j = 2 we need the binary expansion of 2^ — 2 = 0 which 
has zero ones. Thus we get that = —1 and = -|-1. 


(ii) We assume the lemma holds for n = k, i.e., for row k of the tree graph in Fig. Hence, for 
J = l,2 ,..., 2^ 

fi,k) _ ^_2^^Par{/3o,...,/3fc-i}^ 

From the tree presented above we observe that the 2” signs of an arbitrary row are the same as 
the last 2" signs of row n + 1. 

Consider now the (fc + l)-th row of the tree. Using the induction h5rpothesis and the observation 
above it follows that the lemma holds for the last 2^ signs of the (k + l)-th row as well. We can 
also see this by observing that for the last 2^ signs of the {k + l)-th row we are interested in 
the binary expansions of 2^+^ — j for j = 2* + 1 ,..., 2^+^, which is essentially equivalent to 
considering the binary expansions of 2^ — j for j = 1,..., 2^. This shows that, for j = 1,..., 2^, 

= ^a+2hfc+i)^ 


What remains is to prove the lemma for the elements 2^ — j, j = 1 ,..., 2^, of the (fc + l)-th row. 
At this point, we observe that at an arbitrary row n, because of symmetry 

_ _c0+2""h")_ 

Hence, the signs of terms j and j + 2”“^ in the n-th row will always be opposite. This yields 
that in the {k + l)-th row we have, for j = 1 , 2 ,..., 2 ^, 

^O'.fc+i) ^ _^(j+2hfe+i)^ 
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But we know that 


g(j+2'‘,fc+l) _ ^•_2^^Par{/3o,..../3fc_l,0} 


We also know that, for j = 1,2,..., 2^, the binary representation of 2^+^ — j has one more 1 fhan 
fhe binary represenfafion of 2^+^ — 2^ — j = 2^ — j. This leads fo fhe expression 


and in addifion 


for all j = 1, 2 ,..., 2'=+h 


(_l)PaT{/3o,...,/3fc_i.l} _ _j-_2^Par{/3o,...,/Sfc_i,0} 


□ 


Before proceeding fo fhe proof of Thm. 3.2 we make some remarks concerning fhe resulf esfablished 
in Lemma 3.2 fhese remarks are used in fhe proof of Thm.|3^ 

Remark 4. For an arbifrary row n in fhe free presenfed in Fig.|^ we have fhaf 

^(1,") = _c(h"+i)^ 


We know fhaf in order fo find fhe firsf sign of fhe n-fh row we musf find fhe binary expansion of fhe 
elemenf 2” — 1, which has exacfly n ones. Thus, for an arbifrary n > 1, we gef fhe expression 

cihn) ^ 

and fhis also shows fhe relafion menfioned in fhe remark. 

Remark 5. For I = 1,..., n and j = 1,, 2"“* we have 

^(2h-2'-i + l,n) ^ _^(2h-2'-i + l.n+l)^ 


To see fhis we observe, as in fhe proof of Lemma 3.2 fhaf fhe 2" signs of fhe n-fh row are fhe same as 
fhe lasf 2" signs of fhe (n + l)-fh row. This gives, for I = 1,..., n and j = 1,..., 2"“*, 


^(2h-2'-Wl.n) ^ ^(2h-2'-W2"-Wl.n+l)^ 


(5.8) 


Using fhe symmefry of fhe signs in each row, i.e., for i = 1,..., 2" we obfain 

fhe equalify above. 

Remark 6. For fhe j-fh sign of fhe n-fh row and fhe (2" + j )-fh sign of fhe (n -I- l)-fh row we have fhe 
following expression, 

.. (5.9) 


_ ^(2"-|-j,n+l)_ 


For fhe value of we need fhe binary represenfafion of 2" — j while for fhe value of + 3 ,n+i) 
we need fhe binary represenfafion of 2”+^ — 2 ” — j = 2” — j, giving 15.9 1. 
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Proof of Theorem ^?2\ Now, having Lemma 3.2 at our disposal, we proceed with the proof of Thm. 
|3.2| Consider fhe case n = 1. Then we have one exponentially distributed random variable Ti with 
parameter qi. Thm. 3.2 gives 

P.(Qri Gdy) = + 


From Lemma 


3.2 


we have that = —landc^^’^^ = +1. Due to Definition 


3.2 


L 


(2!i) “ +^(9i)e 


-'^in)v 


Hence we obtain 


Px(Qti G dy) = — — y) +'^{qi)e dy. 

This is the expression found in (|2.2|, and we conclude fhaf our resulf holds for n = 1. 


We now assume fhaf Thm. 3.2 holds for n = k — 1. Consider now fhe case oi n = k exponenfially 
disfribufed random variables, fhen conditioning on fhe value of fhe workload af fhe firsf k — 1 expo¬ 
nential epochs yields 


Px(QTi+...+Tfc G dy) 

pOO 

= / P^(QTi+...+rfc_i G dz)¥^{QT^ e dy) 

Jz^O 

= j (x-z) 

•^^=0 \ i=l 


fe_l 2'"-*-! 


+ E E Gdy)d 

1=1 j=l ) 

k1 pCXi 

hi.fc-i) Y[qi * ... * {x - z) (-qkW^‘^'‘\z - y)) dz 

i=l "'^=0 

qi-qj(q^^)e-'^(<i'^)v * ... * (x - z)Z^’i>‘\z)dz 

O 

■Z-(2E^2'-1+I.i)(^) {^-qkW^‘^>’\z - y)) dz 

* ... * (x) 

We see fhat we have fo evaluafe fhe following four infegrals: 

Ji = ^ (x-z)W(«'')(z-y)dz, 

Ja = ^ (x-z)z(«'=)(z)dz, 


2=1 

k — 1 

-^E E 

1=1 j=i 


I z=0 

J-i = —9fe 


/2=0 


-(fc-i) 

■^(2h'-2*-i+1.0 


(z)W^‘^''^(z — y)dz, 


/ OO 


dy. (5.10) 

(5.11) 

(5.12) 

(5.13) 

(5.14) 
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These integrals J7i,..., J 4 can be interpreted by looking at the tree graph in Fig.j^ J\ is related to the 
scenario we are at node 1 of row k—1 and we move down and left, while J 2 is related to the scenario 
we move down and right. suggests that we are at some node in row fc — 1 which lies in a subtree 
with root k — l — 1 rows above the node under consideration and we move left, while corresponds 
to the case we move down and right. 


o Integral J\. By a change of variable argumenf and using fhe facf fhaf = 0 for a: < 0, we 

find fhaf 

(5.15) 


* ... * ( 3 . _ yy 


From 15.10 1 we see fhat fhe firsf ferm of f'a:{QTi+...+Tk € dy) is equal to 

k 

* ^( 91 )^ (x-y). 

Using Remark|^we have and this shows that we have identified fhe firsf ferm of 


the expression in Thm. 3.2 


o Integral J 2 . It is straightforward fhaf 

J 2 = * ... * (x). 


(5.16) 


Hence, fhe second ferm of 15.101 is equal fo 


fc-i 


..(i.fe 1 ) (?i'F(gfc) exp[—T'(gfc)y] * ]4/(9fe-i) * ... * (x). 


This ferm corresponds fo Z = fc in fhe summafion of fhe second ferm in Thm. 3.2 We need fo show 
fhat 


. (fc) 

^(2'=-i + l,fc) 


fc-1 


= Y[ (Z,«'(qfc)exp[-T'(gfc)y]- 


Z=1 


First we observe that, for I — k and j = 1, we have fhat to( 1, fc) = fc and due to Remark]^ we have 
that This shows that 


k — 1 

di.fc-i) g^^(^^) 
2=1 






2=1, 

27 ^m(l,A:+l) 


yielding that the expression for fhe second ferm, for I = fc, agrees wifh Thm. 3.2 


The infegrals jTa and J 74 correspond fo convolufion ferms in one of fhe subfrees mentioned before 
(see Fig.0. Therefore, fhe order of fhe convolufion ferm does nof change and fhese infegrafions only 
change fhe coefficienfs L. 

o Integral J 3 . We now show fhaf equals t| 2 **^_ 2 !-i+i /)(y)- Fig.j^we see fhaf when we move fo 
the left the numbering of fhe coefficienfs remains fhe same. When looking in row fc — 1 at the element 
^( 2 q- 2 '-i+i i)' ^ “ 1 ,..., fc — 1 and j = 1,, 2 ^“*“^, and moving down and left we arrive at the 

(k) 

element L]^ 2 ‘j- 2 ‘-^+i i) which has the same 'labeling'. 
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Since we have an expression for h| 2 'j- 2 '-i+i i) induction hypothesis yields 


Js = -c^ 


k rn.{j,l) ^ k-2 


n 


J.JL ±± q^-Qi+i, 

2=1, I—I 2=m(j,/) + l ’ 


f'OQ 

/ — j/)dz. 

Jz=0 


(5.17) 


By a change of variable and the definition of the g-scale function (•) in Eqn. (2.1 1 , it follows that 

1 

f z—0 Q7n{j,l) Qk 


pOO 

J z^O 


)y. 


(5.18) 


Substituting Eqn. ( 5.18| into Eqn. ( 5.17| , in combination with the use of Remark|^in the second step, 
we find 


k rn{j,l) 






i=l, i=l i=mU,l) + l 




(5.19) 


as desired. 


Integral J 4 . Before proceeding to the last integral, it is noted that from Eig. we observe that 

, for some I = 1,..., k — 1, j = 1 ,..., and row k — 1, and 


when we are at node h| 2 q- 2 '-i+i i) 


we move down and right, then we obtain the coefficient ly But this is equivalent to 

^( 2 q- 2 '-i+i i) i ~ + 1,..., 2^~K Erom the definition of the Z-scale function in Eqn. (2.2 


and interchanging integrals, it follows that 

pOO 


1 


qm{j,l) 


4'(g 

771 qm(j,l) Qk 

Due to the fact that j = 1,..., and / = 1,..., fc — 1 we have that m{j, 1) < k — 1, this leads to 

fe-i 


Ja = c 




1 


k-1 m(j,l) 

n® It n 

L, L, «-®+'i.4h+i 


(5.20) 


The last property we should verify is that, for I — 1,... ,k — 1 and j = 1,...,2* ^ it holds that 
J 4 = h| 2 i*j_ 2 i_i^^^ 2 '"-i /)■ But we observe that this is equivalent to showing that, for I = I,... ,k — I 
andj = 2'=-'-i + l,...,2'=-' 


J4 = L 


(k) 


(5.21) 


^(2q-2'-i+i,i)' 

Eor the values of I and j we consider, we have that m{j, 1) = k and thus, by Remark]^ we see that 
Eqn. ( 5.20| gives Eqn. ( 5.21| . 
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Returning to <5.10|, and using the expressions found in |5.15|, (5.16\, (5.19| and <5.21|, we find 


lPa:(QTi + ...+Tfc G dj/) 
k 

= Y[ q, * ... * (a; _ y) 

i^l 

k 

Y[ *... * ( 3 .) 

2=1 


k -1 2 '® 




1=1 j=l 


dy. 


The lasf expression can be wriffen more compacfly, yielding 


^x{QTi+...+Tk G dy) — 


K, 

]J qi * ... * (a; _ y) 

2 = 1 

+E E ^(2E2 -^+i.o (y) * • ■ ■ * 


1 = 13=1 


dy. 


This concludes fhe proof of Thm. 3.2 


□ 


Proof of Corollary ^l\ The expression in Corollary |3.1| is derived as a sfraighfforward application of 
fhe resulf esfablished in Thm. 3.2 For fhe friple fransform we know fhaf 

poo poo poo 

/ e-^“E^e-“®"'i+'■+^"da; = / 6"^“ / 6"“^P;,(Qri+...+T„ G d2/)dx. 

Jo J x—0 J y—0 

Using Thm. 3.2 we see fhaf 

poo poo poo 

/ e-'^^E^e-“®^i+ +^"dx = c(i>")TTg, / / e-^^e-°‘y k {x - y)dydx 

Jo Jx^o Jy^O ^ ' 

/ \ Z*'^ 


-^EE 

/=i i=i 


' x—O 


(n) 


where fhe coefficienfs T( 2 'i- 2 '-i+i i) 
following fhree infegrals: 


are given in Definifion 


3.2 


We see fhaf we have fo work wifh fhe 


/Cl = 


00 pOO 


a-p^^-o^y 


* ... * [x - y)dydx, 

IC 2 = dx, 

00 


X—0 'J y—0 
00 


1C. = 


-ay T {n) 
O k'/.l 


y=0 


(2<a-2'-i + l,0 


{y)‘iy- 


Relying on fhe properties of fhe g-scale funcfions and affer some sfraighfforward calculus, we find 


/Cl = 


n 


a + P $(/3) - q^ ’ 




P fj{ HP) - 
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and 


JQ = c(2‘,-2—+ l,n 




T-U,l) 

n n 


n 


^ ^{Qm(j,l)) ■ -I ■ / / ■ . , Qi Qi-\-l - / ■ I 1 ^*+1 


We conclude that by substitution of /Ci, IC 2 , and /C 3 in < 5.22 1 , we have established Cor. 3.1 


□ 


6 Conclusion and Discussion 


In this paper we have analyzed the transient behavior of spectrally one-sided Levy-driven queues. 
We considered the joint behavior of Qti, Qti+Ts, ■ • ■, Qti+...+t„ where Ti is exponenfially disfribufed 
with parameter qi, and we specifically focused on From fhe main resulfs if follows fhaf 

fhis fransienf behavior obeys an eleganf and appealing free sfrucfure. Inferesfmgly some numerical 
illusfrafions showed fhaf E^; Qt is firsf decreasing in t and fhen converges fo fhe sfeady-sfafe workload 
from below in case x is chosen 'slighfly' above fhe sfafionary workload. 

We have resfricfed ourselves fo analyzing Qt wifh T disfribufed as fhe sum of n independenf expo¬ 
nential random variables, but our result is readily extended to that of Qt wifh T obeying a Coxian 
disfribufion. This is a parficularly useful facf, as any disfribufion on fhe posifive half line can be ap- 
proximafed arbifrarily closely by a sequence of Coxian disfribufions, see e.g. 121 Section III.4]. In more 
defail, fhe analysis looks as follows. Consider fhe sifuation fhaf T follows a Coxian disfribufion wifh 
n phases; we lef fhe lengfh of phase i be drawn from an exponenfial disfribufion wifh paramefer qi, 
and we lef fhe probabilify of moving from phase i io i + 1 he pi (wifh fhe convenfion fhaf = 0 ). 
Then, for fhe specfrally-posifive case. 


fc-i 


E,^ 


„ —a Qt 


= -Ffc) 


^-a QT-^ + ...+Tf. 


k=l 


3.1 


The densify in fhe specfrally-negative case follows 


where E^; ®^i+ is as obfained in Thm. 

by a similar argumenf. 

To conclude, we like fo menfion some topics fhaf are of inferesf for fufure invesfigafion. Alfhough 
fhe class of Coxian disfribufions for fhe epoch T is sufficiently rich, it might of inferesf fo sfudy 
fhe behavior of Qt if T has a general phase-fype disfribufion. Specifically, we did nof explicifly 
derive fhe resulfs in case some parameters qi are identical. This follows as a direcf application of 
THopifaTs rule, buf fhe expressions fend fo become cumbersome. Anofher open question concerns 
fhe fransienf behavior for specfrally fwo-sided Levy processes. Finally, we expecf fhaf fhe fransienf 
analysis presented here may be applicable in inference procedures, fo esfimate fhe queue's Levy inpuf 
process from a finife number of successive workload observations. 
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